Homology 


Because the physics which goes into determining the structure of main- 
sequence stars does not change rapidly with mass, one can make the 
assumption that a star with mass Mı will just be a scaled version of 
a star with mass Mo. Stars built under this assumption are called ho- 
mologous. Through homology, relations can be found describing how 
various stellar quantities (such as radius, luminosity, central tempera- 
ture, etc.) change with mass. 


To begin with, consider that if star 1 is really just a scaled up version 
of star 0, then the radius and mass at any point in star 1 will be 
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where the subscript 0 refers to the reference star. From the above 
relations, it is obvious that 
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With this relation, we can compute the density at any point in the 
star. From the equation of mass conservation, 
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Similarly, we can use the equation of momentum conservation and the 
assumption of hydrostatic equilibrium to compute a relation between 
pressure, mass, and radius 
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The trivial integration then gives 
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To make further progress with homology, let’s assume that energy gen- 
eration, opacity, and the equation of state are all power laws of density, 
pressure, and temperature, and that the power law exponents are the 
same for every star. Under this assumption, and the definition of the 


thermodynamic quantities (1.5), (1.6), and (1.7), the stellar densities 
will be given by 
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Using (11.1) and (11.14), the nuclear reaction rates will be 
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and, via (9.13), the opacity can be parameterized by 
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In practice, of course, these exponents will not be the same in every 
star. For instance, the temperature dependence of nuclear energy gen- 
eration in a low mass star is significantly different than in a high-mass 
star. However, for comparisons of similar mass stars, this approxima- 
tion should be reasonable. 








Now consider the equations of mass conservation (2.1.3), pressure bal- 
ance (2.2.4), energy conservation (2.3.2), and thermal structure under 
radiative energy transport (2.4.4) and (3.1.6). From (17.3), mass con- 


servation is 
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which, through the equation of state (17.5), becomes 
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From (17.4), pressure equilibrium is simply 
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The relation of energy conservation is also straightforward 
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If we perform the trivial integration, and substitute (17.6) in for €, 
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From the equation of state (17.5), we then get 
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Finally, we can write down the equation for the thermal structure of 
the (radiative) star 
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When we substitute for r/rp using (17.1), and integrate over dT and 
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If we substitute for the opacity using (17.7), then 
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Equations (17.8) - (17.11) now allow us to examine how the radius, 
pressure, temperature, and luminosity of homologous stars change with 
mass and chemical composition. First, to study the behavior with 
mass, set u = uo. Equations (17.8) - (17.11) in matrix form are then 
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Solving these four equations with four unknowns is tedious; the result 
is 
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Similarly, by holding the mass constant and allowing p to vary (i.e., 
by replacing the right-hand column vector of (17.12) with the vector 
[—y, 0, Ay, 0|), one can solve for the dependence of stellar structure on 
composition: 
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Note that the solutions to the homology equations are complicated. 
In general, stellar structure depends on the nuclear reaction rates and 
opacities and the equation of state. However, by substituting in for 
reasonable values of a, ô, A, v, s, and t, one can get an idea of how stellar 
properties should scale with small changes of mass and composition. 
For example, if 6 = 0 (i.e., a polytrope with n = a/(1 — a)), we have 
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which is the identical expression derived for the mass-radius relation 
of a polytrope. 


HOMOLOGY RELATIONS WITH THE IDEAL GAS LAW 


Because homology relations are derived using ratios, the pro- 
cedure cannot be used if any of the terms are additive. Thus, a 
general equation of state which include gas and radiation pres- 
sure cannot be modeled in this manner. However, if we restrict 
ourselves to ideal gases, and use reasonable assumptions for the 
energy generation rate and opacity, then we can use homology 
to predict stellar structure. 


For ideal gases, a = 6 = y = 1. Moreover, most nuclear reac- 
tions involve two-body collisions, hence the energy generation 
per unit volume goes as p°. Thus, à, the density coefficient per 
unit mass for nuclear reactions, should also be A = 1. 








We can now investigate the homology relations for the three 
opacity laws (electron scattering, Kramers, and H7), and for 
values of v between ~ 5 (for the proton-proton chain) and ~ 20 
(for the CNO cycle). The results are shown on the next page. 





Dependence of Stellar Structure on Mass 


Opacity P=o-. he TY) 


log R — log M Coefficient 
e-scattering 0.50 0.69 
Kramers 0.20 0.52 
H- 0.72 0.78 


log P — log M Coefficient 


e-scattering 0.00 —0.77 
Kramers 1.20 —0.08 
H- —0.89 —1.13 


log T — log M Coefficient 
e-scattering 0.50 0.31 
Kramers 0.80 0.48 
H- 0.28 0.22 


log £ — log M Coefficient 
e-scattering 3.00 3.00 
Kramers 5.40 5.24 
H- 122 1.83 


log Teg — log M Coefficient 
e-scattering 0.50 0.40 
Kramers 1.25 1.05 
H- —0.92 0.05 


log p — log M Coefficient 


e-scattering —0.50  —1.08 
Kramers 0.40 —0.56 
H- —1.20  —1.36 
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0.34 
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—0.98 
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A few items to note: 


e For electron scattering, the mass-luminosity relation for stars 
is independent of the mode of energy generation. The relation 
is entirely defined by hydrostatic equilibrium, the ideal gas law, 
and opacity. Thus, stars in this regime will adjust themselves 
so that the nuclear reactions will provide the requisite amount 
of luminosity. 





e The mass-radius relation for all relevant coefficients is a small 
positive number. Thus, unlike n = 3 polytropic stars (which 
have a coefficient of —1/3), the radius of homologous stars will 
increase (slightly) with mass. 





e Since Teg is formed simply from £ and R, the dependence of 
Tex On mass is easily found. The coefficient is positive, but the 
slope of the relation depends on the opacity source. Similarly, 
we can derive the slope of the relation in the log £ — log Teg 
(HR) diagram. 


Opacity v=o. Pali vals w= 20 


log £ — log Teg Coefficient 
e-scattering 6.00 7.43 8.31 8.90 
Kramers 4.32 4.99 5.36 5.60 


e In a similar manner, we can calculate the behavior of £ and 
Teg with mean molecular weight, u. Once again, in the electron 
scattering case, the dependence of £L on p is independent of 
nuclear energy generation. In fact, (17.14) gives £L « ut; the 
luminosity is more sensitive to u than it is to mass! 





Dependence of Stellar Structure on Molecular Weight 





Opacity bes yal ysl v = 20 
log R — log u Coefficient 

e-scattering ee 0.46 0.61 0.70 

Kramers —0.47 0.12 0.37 0.51 

H- 0.50 0.61 0.68 0.73 
log P — log u Coefficient 

e-scattering aie —1.85 —2.44 —2.78 

Kramers 1.87 —0.48  —1.49 —2.04 

H- —2.00 —2.43  —2.71 —2.91 
log T — log u Coefficient 

e-scattering — 4.23 5.28 5.87 

Kramers —2.27 1.84 3.60 4.58 

H- 4.50 5.26 5.75 6.09 
log £ — log u Coefficient 

e-scattering Ja 4.00 4.00 4.00 

Kramers 8.73 8.44 8.31 8.24 

H- 1.00 2.09 2.79 3.27 
log Tog — log u Coefficient 

e-scattering Pets 0.77 0.69 0.65 

Kramers 2.42 2.05 1.89 1.81 

H- 0.00 0.22 0.36 0.45 





e We can also calculate the dependence of Tyg on u, and there- 
fore the resulting u vector in the HR diagram. The resulting 
slope is shallower than the main sequence, but greater than the 


lines of constant radii. 


Opacity v=%5 g= 


log £L — log Teg Coefficient with u 
e-scattering 4.27 5.20 
Kramers 3.43 4.04 


5.76 
4.35 


6.13 
4.53 


